Efficient numerical method to calculate three-tangle of mixed states 
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We demonstrate an efficient numerical method to calculate three-tangle of general mixed states. 
We construct a "energy function" (target function) for the three-tangle of the mixed state under 
certain constrains. The "energy function" (target function) is then optimized via a rephca exchange 
Monte Carlo method. We have extensively tested the method for the examples with known analytical 
results, showing remarkable agreement. The method can be applied to other optimization problems 
in quantum information theory. 
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One of challenges in present quantum information the- 
ory is the quantification of quantum entanglement. The 
primary studies on entanglement measures focus on bi- 
partite systems. In bipartite entanglement measures, for- 
mation of entanglement, presented by Bennett et al., is 
very fundamental [l!] . The significance of this measure re- 
lies on not only be able to provide an accurate boundary 
between separable states and entangled states but also 
constructing a general formula for entanglement mea- 
sures of mixed states: once a measure for pure states 
is given, the corresponding measure for mixed states can 
be obtained via the convex-roof extension. This scenario 
can be even generalized to quantify multipartite entangle- 
ment. The initial attempt is 3-tangle of 3-qubit system, 
provided by Coffman et al Entanglement measures 
beyond 3-qubit were also presented [1, |j] . Although such 
kind of definition is straightforward, the practical eval- 
uation of the convex-roof extension is a difficult mathe- 
matical problem. To date a general analytical method is 
known only for the concurrence of two-qubit mixed states 
[5, 6]. 

Concurrence of mixed two qubits has been applied to 
study quantum phase transitions. Most of studies show 
that, for general spin 1/2 lattice models, pairwise entan- 
glement (or its functions) depicted by concurrence of the 
nearest-neighbor two site has special singularity at quan- 
tum critical points (see[2j and its references). A natu- 
ral generalization for this viewpoint is that multipartite 
entanglement beyond two sites will reveal more exten- 
sive characteristics on quantum phase transition. How- 
ever, as far as 3-tangle is concerned, since the convex 
roof is obtained by minimizing the average tangle of a 
given mixed state over all possible decompositions of pure 
states, there is no general method to calculate the three- 
tangle of a mixed state. The analytical method is only 
capable to study some particular three qubit states |8l-tlll|. 

In this work, we develop a method to calculate the 
three-tangle of any mixed states by numerically mini- 
mize certain target functions. The target functions are 



ogtimized by using replica exchange Monte Carlo (MC) 
[H, [l^l method. Replica exchange MC has been widely 
used in condensed matter physics to simulate the systems 
with very rough "energy surfaces", such as spin glasses 
[T3 | , protein systems and Lennard- Jones particles [11, [ll| 
etc. The replica exchange MC method is much more ef- 
fective to find the global minium of the target functions 
than the conventional simulating annealing method. By 
applying this method we have successfully obtained the 
three-tangle of mixed (generalized) GHZ and W states. 
We also obtained the three-tangle of GHZ state with 
white noise, whose density matrices are of rank 8. 

The concurrence C{iPab) measures the biparticle en- 
tanglement between particle A and B in a pure two-qubit 
state \iPab} ■ C is defined as 



(1) 



(f>ii}of \iPab) 
The measure 



in terms of the coefficients {(t>oo, (j^oii 4'io 
with respect to an orthonormal basis, 
for three particle entanglement in a three-qubit state 
\ipABc} has been introduced in Ref.Q as three-tangle 
Tsi^ABc)- It can be expressed in terms of the coeffi- 
cients {0000, ^^'001, 011l}- 
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d2 = 0OOO01110O11</'1OQ + 0OOO0111</'1O10O1O 
+0OOO0111011O</'OO1 + 0O1101OO</'1O10O1O 
+0O1101OO011O</'OO1 + 01O10O1O</'11O0OO1 (4) 

^3 = 0OOO</'llO</ilOl'/'Oll + '/'1110OO1</'O1O'/>1OO (5) 

It can be shown that for any factorized state, the three- 
tangle vanishes. For the GHZ state. 
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Tj,{GH Z)=l. It also appears a class of entangled three- 
qubit states represented by \W) state with vanished, 



I.e. 



\w) = ^(|ioo> 



|010) + |001)) 



(7) 



The three-tangle of mixed state can be obtained via 
convex-roof extension jll. Isl. fl7l| . Suppose the density ma- 
trix of a mixed state p can be decomposed into sum of 
some pure states tTj, i.e., 



where, 



(8) 



(9) 



is the wave function of a normalized pure state. The 
three-tangle of the mixed state is defined as the average 
pure-state concurrence minimized over all possible de- 
compositions. 



r3(p) = min^p^T3(7ri 



(10) 



It is very difficult to develop an universal analytical 
method to calculate the three-tangle of mixed three- 
qubit state. So far, there are very limited examples 
of mixed states whose three-tangle have been obtained 
analytically [l-[ll|. Alternatively, one could resort to the 
numerical methods. To get the three-tangle of mixed 
states is a constrained minimization problem. The pure 
state wavefunction can be expanded on the orthonormal 
three-particle basis, i.e., = Y^^Cia\^a)- The first 
constrain is, 



Pi Ci,QCjj3 — Pa/3- 



(11) 



where is the density matrix in the basis of {j^a)}. 
Furthermore, the coefficients have to satisfy the normal- 
ization conditions. 



K is a large constant, to ensure ^Q. iVc = 8 is the 
dimension of three-particle basis set. Np is the number 
of partition in the decomposition. We can easily see that 
Np determines the number of variables = Np+2NcNp. 

In EgfTHl K should be large enough to ensure a small 
residual to numerically satisfy the constrain given in 
Eq llll However, large k will make the "energy" surface 
of E very rough, i.e., has many local min- 

ima separated by large "energy" barriers, which causes 
great difficulties in finding the global minima of three 
tangle using traditional constrained simulated annealing 
(CSA) method, because it tends to be trapped at some 
local minima. Here, we adopt the replica exchange (also 
known as parallel tempering [ll|) Monte Carlo method 
fl2| which simulates M replicas simultaneously each at 
a different temperature /3o = l/Tlnax < /3i < ... < 
l3n~2 < Pm-1 = l/7niin Covering a range of interest. 
Each replica runs independently, except that after cer- 
tain steps, the configurations can be exchanged between 
neighboring temperatures, according to the Metropolis 
criterion, 



w 



1 AH <0 
e^^^ otherwise 



(15) 



where AiJ = —{jii — j3i-i){Ei — i?i_i), in which Ei and 
Ei_i are the energy of the i-th and i-l-th replica. During 
the exchange, the detailed balance in canonical ensemble 
is satisfied. Importantly, the inclusion of high-T config- 
urations ensures that the lower-T systems can access a 
broad phase space and avoid becoming trapped at local 
minima. The replica temperatures are adjust so that the 
exchange rate between the replicas are all about 20% p^. 
while keep the highest temperature /3o = l/Tlnax and 
lowest temperature /3m-i — ^/T^i^ fixed. This can be 
done by adaptively adjusting the temperature of each 
replica using a recursion method proposed by Berg[20l|. 
We slightly revised the the algorithms to give better per- 
formance. Suppose, in the n-th iteration, the temper- 
ature of the i-th replica is (if and the acceptance rate 
between /3" and is a['". In the n-|-l-th iteration. 



^ \cia\^ = 1, and = 1,Pj > . (12) /Jf"^^ can be updated as. 



The above constrains can be enforced in the simulation 
by explicitly applying normalization factors at each MC 
step, whereas the constrain Eq. ^2 can be enforced via 
a penalty function. The "energy function" (target func- 
tion) reads. 



where, 



n+l 



/3,r-Y + (i-c+car)(/3r-/3r-i), 



M-l 



(16) 



(17) 



^(fe,^',}) = ^p,r3(*.) + Ki?2 



(13) 



where {p^, realizing the mixed state as given in Eql^l 
is the residual between the searched density matrix 
and the target density matrix p, defined as 

i=l 



c is an empirical parameter that can be adjusted to ac- 
celerate the convergency. In the present case, we found 
that by choosing proper c < 1 (in this case c=0.7) the 
recursion scheme get a suitable /3 distribution within tens 
of iterations. We then fixed the replica temperatures for 
the simulations. 

Based on convex-roof extension, Lohmaycr et al. have 
provided a complete analysis of mixed three-qubit states 
composed of a GHZ state and a W state, obtaining 



the optimal decomposition and convex-roof for three- 
tangle[8j. Eltschka et al generalized this method to treat 
the three-tangle of mixed three-qubit states composed of 
a generalized GHZ state and a generalized W state. The 
explicit expressions for mixed-state three- tangle and the 
corresponding optimal decornposition for the more gen- 
eral form were also presented[3|. We first test our scheme 
by comparing with the available analytical results. We 
construct a series of target functions based on the density 
matrix given in Refs. as were described in the previ- 
ous sections. Replica exchange MC is then performed to 
find the three-tangle and optimal decomposition which 
are compared with the analytical results. 

In all our minimization process, the temperature range 
is chosen between T-a^^^=lQQ and T,„in=10~^, such that 
Tmax is high enough to help the low temperature replica 
to avoid being trapped in the local minimum and Tmin is 
low enough to give the final minimized value an accuracy 
of 10^^ theoretically. With fixed temperature range, the 
suitable number of replicas is determined by the corre- 
sponding density matrix and weight factor k. Usually 150 
replicas are enough. In our simulations, k is adjusted to 
be lO"* - 10"^ to keep the final < 10-^° which is con- 
sidered accurate enough in identifying the three-tangle of 
mixed state. In each simulation, the number of partitions 
Np is fixed. We then increase Np until the calculated 
three-tangle converges. 

Let us first look at the mixed GHZ and W states, The 
density matrices of the states are, 

p{p) = pTTGHZ + - p)'^w , (18) 

where ttghz = \GHZ){GHZ\ and tt^ = \ W){W\ are the 
density matrices of GHZ and W states respectively [8]. It 
is known that t^{GH Z)^\ and T3(W^)=0. According to 
Caratheodory's theorem, Np=A is enough to get the op- 
timized Ts of a mixed states of rank 2Q. We calculate 
by using Np=A - 8. Indeed, we find in our simulations, 
that Np=A is enough to give the optimal values. The re- 
sults are shown FigUJa), which depicts the numerically 
calculated T3 of the mixed GHZ and W states (open cir- 
cles), compared to the analytical values (solid line) 8]. As 
one can see that the numerical and analytical results are 
in excellent agreement. 

We then calculate r3 for the mixtures of generalized 
GHZ and generalized W states whose density matrices 
are, 

p{p) = p \gGHZ) {gGHZ\ + {1 ~ p) \gW) {gW\ , (19) 

where the density matrix of a generalized GHZ state is 

\gGHZ) = a|000) -I- 6|111) , (20) 

and density matrix of a generalized W state is 

\gW) ^ c\00l) + d\010) + f\100) . (21) 

The coefficients must satisfy the normalization condition 
lap + |6p=l, and \c\^ + jdp + \f\^=l. For the pure gen- 
eralized GHZ and W states, it is easy to calculate that 
T3{gW)=0, and TsigGHZ) = 4\a^^\. 
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FIG. 1: The three-tangle calculated with replica exchange 
compared to analytical results, for (a) the mixed GHZ and 
W states, (b) the mixed generalized GHZ and generalized W 
states, (c) the mixed GHZ state, W state and flipped W states 
with n—2, (d) the mixed GHZ state and white noise states. 



TABLE I: Numerically calculated T3 of mixed GHZ, W and 
flipped-W states for n=2 compared to the analytical results 
given in Ref . [lol . 



p 0.2 0.75 0.8 0.85 0.9 0.933 0.95 

analytical 0.1835 0.3775 0.5805 0.7182 0.7897 
numerical 5.7e-6 2e-5 0.1836 0.3781 0.5811 0.7186 0.7901 



The calculated rs of the mixed states are shown in 
FiglUb) for a=0.2, c=0.2, d=0.2. As we see the results 
are in excellent agreement with the analytical results. 
Again we find that Np=A is enough to the optimal results. 

In the above tests, we demonstrate that our method 
is very effective to deal with the three-tangle of rank-2 
mixed states. Recently, Jung at al. has obtained the 
analytical results of the three-tan gle o f the mixed states 
of GHZ, W, and flipped-W states[10|. These states are 
of rank 3. The flipped-W states is defined as 

|u;) = -^(|110) + |101) + |011)). (22) 

The density matrix of mixture of GHZ, W, and flipped-W 
states are given by, 

Pip, q)^p \GHZ) {GHZ\+q \W){W\ + il-p-q) \W) {W\ , 

(23) 

where q was defined as q — 1—p/n, n is a positive number. 
We then compared our numerical results to the analyti- 
cal values. The results for the case of n=2 are shown in 
FigHJc). The numerical results are also in good agree- 
ment with the analytical results. We further calculate 
three-tangle of the mixed states which are of rank 4, given 
in Ref|Tl|, obtaining the same accuracy compared to the 
analytical results. 

To give an illustration of the accuracy, the numeri- 
cal values of mixture of GHZ, W and fiipped-W state 
with n=2 are also listed in Table U We can see that 
all the results are within an accuracy of ~ 10~^. The 
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FIG. 2: (Color online) The A'p dependent Ara calculated with 
replica exchange for density matrix constructed by GHZ state 
with white noise, where Ars = T3 - r™'" with r™'" the min- 
imum three-tangle for each p value. 



accuracy can easily be improved by increasing k and de- 
creasing Tniin. For the example given above, in the case 
of p=0.8, if we reduce Tj^in to 10"-^° and use k = lO^*^, we 
obtain T3=0. 18349906 compared to the analytical value 
0.18349861 in an accuracy ^ 10^^. In principle, our nu- 
merical method can achieve arbitrary accuracy. We have 
also tried to optimize Eq. [T3]via traditional CSA method. 
We find the performance of the traditional method is very 
bad, because it can easily be trapped at some local min- 
ima. 

Having demonstrated the ability of the method, we 
then calculate the T3 of the GHZ states mixed with white 
noise, which does not have analytical solutions so far. 
The density matrices of the states are 



p{p) = p TTGHZ + T^E , 



(24) 



where tt^; is a unit matrix. In this case, the density matri- 
ces are of rank 8. We need at least Nr,=^ to ensure ~ 



0. In practice, we have calculated up to Np=2Q. We 
find that use 7Vp=15 can converge to less than 10~^ 
for this problem. (See Fig. I^). The calculated results 
are shown in Figlljd). From Fig. [TJd), we see that 
of the white noise mixed GHZ states is zero for p less 
than about 0.7, then increase almost linearly to 1.0 as 
p increases. The behavior of looks similar to that of 
mixed GHZ and W states. 

It is straight forward to apply the method to any mixed 
states. It therefore opens up a way to study the multi- 
particle entanglement effects in many important systems, 
e.g., in systems with quantum phase transitions, which 
was impossible before, because the analytical solution to 
the tree-tangle of a general mixed state is not available. 

To conclude, we have demonstrate an efficient numeri- 
cal method to calculate the three-tangle of general mixed 
states. We construct a "energy function" (target func- 
tion) for the three-tangle of the mixed state under cer- 
tain constrains. The "energy function" (target function) 
is then optimized via a replica exchange Monte Carlo 
method. We have tested the method for the examples 
with known analytical results, showing remarkable agree- 
ment. We further calculate the three-tangle of GHZ 
states mixed with white noise. The work opens a way 
to study three-tangle of general mixed states, and can be 
generalized to solve many other optimization problems 
in quantum information theory. 
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